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Abstract 

Interrelations of domain wall (extrinsic) contributions to dielectric, piezoelectric and mechani- 
cal properties of ferroic samples are a hot subject of both theoretical and experimental research. 
Recently, we have derived theoretical formulas for such contributions, using model of composite 
layer sample: central single-crystal ferroelectric-ferroelastic layer, isolated from electrodes by pas- 
sive layers. Here we present more general results, discussing the dependence of the contributions 
on geometric and material parameters of the composite and including in a special case intrinsic 
piezoelectricity. We also discuss the above-mentioned interrelation and compare our results with 
measurements of all the contributions on the same RDP sample in a wide temperature interval 
under the phase transition. It seems that both calculated and observed results remind of the so 
called Pippard-Janovec thermodynamic relations. 
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I. INTRODUCTION 



The problem of extrinsic contributions to dielectric, piezoelectric and mechanical proper- 
ties has been addressed by many authors. In prevailing number, only extrinsic permittivity 
has been studied. Arlt et al. Q] and Herbiet et al. [3] were the first, who addressed 
the domain wall contributions to all involved properties: permittivity e, elastic compliance 
s and piezoelectric coefficient d in piezoelectric ceramics (see also a later work of Zhang 
et al. [4]). In this paper, we have in mind ferroelectric and ferroelastic crystals with only 
two domain states, in particular, the KH2PO4 family. In our recent papers, we studied the 
equilibrium domain structure in thin films [sl and equilibrium contributions to permittivity 



6l, piezoelectric coefficient l3l and to all these properties including elastic compliance 

n p fi 

In these papers |6|, 17|, 18| we have used the model of passive layers (ferroic layer composite), 
theoretically discussed first by Fedosov and Sidorkin , then developed by Tagantsev et al. 



10| and Bratkovsky and Levanyuk 
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12|. 



In this paper we shall concentrate on following questions: how the geometric and mate- 
rial parameters of the composite effect the extrinsic contributions, and what is the mutual 
relation of contributions to e, s and d on the same sample. 

In the following sections, we first introduce the model and basic symbols for geometric and 
material properties of the sample. Then we recapitulate briefly our method of theoretical 
derivation of equilibrium extrinsic contributions, described in detail in jsj. Afterwards, we 
shall summarize and discuss the new results, mentioned above and try to compare them 
with experiment. Also, in a special case of equal material properties of both passive and 
central layers, we include the effect of intrinsic piezoelectricity. These three topics were not 
discussed in |8|. 



II. MODEL DESCRIPTION 



We consider a plate-like electroded sample of infinite area with major surfaces perpendicu- 
lar to the ferroelectric axis z. Domains with antiparallel spontaneous polarization Po,3 = i-Po 
differ in the sign of spontaneous shear eo,i2 = ±eo. We suppose that usual linear state equa- 
tions are valid, including intrinsic piezoelectricity. For simplicity we consider the sample to 
be elastically isotropic. The geometric parameters, characterizing the sample and domain 
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FIG. 1: Geometry of the model in x — z plane. 




FIG. 2: Geometry of the model in x — y plane. 



structure, are shown in Figs. [T]and[2l 

We introduce the following symbols £33^ the component of relative permittivity of the 
passive layer; £33^ the same component for the central ferroic layer; s|^2i2 ^^^1 3^212^ ^^e shear 

component of the compliance tensor of the passive and central ferroic layers, respectively; ^3^2 

(2) 

and d^i2, the piezoelectric coefficient of the passive and central ferroic layers, respectively; 
a+ resp. a_, the width of the corresponding domain; 2a, the period of the domain structure; 
d, the total width of the passive layers; h, the width of the central layer; u, the shear 
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displacement at x = 2a; r = Ti2ext, the external shear stress acting on the sample; V, 
the external voltage on the electrodes; A = (a+ — a_)/(a+ + a_), the asymmetry factor of 
the domain structure, expressing the displacement of the domain walls from the symmetric 
position. 



III. EFFECTIVE PARAMETERS OF THE SAMPLE 

In this short article, it is not possible to reproduce the whole process of derivation of 
effective equilibrium permittivity e'^ at zero external stress, compliance s^fj^g zero voltage 
V, and piezocoefficient df[2 of sample. The process is based on the calculation of the 
total energy density U of the sample, i.e., electrostatic energy, deformation energy and 
domain walls energy per unite surface area of the sample in the x — y plane, as a function 
of free parameters A, V and u resp. r. In equilibrium, the variation of the total energy 
of the isolated system: sample + electric source + mechanical source is zero. Solving this 
standard problem, we can find the equilibrium reaction of the sample on external voltage 
or stress. The process is described in [8| and we shall prepare more detailed publications in 
this direction. 

Because of great complexity of the problem, we shall use the "thick sample" approxima- 
tion, valid when 



/i, (i > a. (1) 

In this case, we can neglect all higher Fourier components of the tensor of deformation and 
electric field, except the constant one. Nevertheless, the formulas for effective parameters 
are still very complex because of the intrinsic piezoelectricity. If this is neglected, we get 
reasonably simple and symmetric results: 



-cff _ A 

'33 ~ ^33 



^) t ("51212'^ + -51212^) + ^0^0^33 d 



^ eg^O (^^33^ + £33/?.^ + Pq (s^i2i2d + ■5m2^ 



(2) 



_ (,) t eleo[ef^d + eSh)+Pis^Si2d 

•^1212 ~ ■^1212 J ~ i^rT^^ „ / /ON TTs T' 



^eff _ fi) Ai)h ep-Ppgo ^ 

"312 — ^^1212 ^33 ^ 2 / (2) (1) \ 2 ( (2) (1) ^ ' ^ ' 
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For thin films, when condition ([T]) does not hold, numerical methods should be applied. 

In the homogeneous composite case £33'' = efr^ = £33, s[2i2 = -51212 — -512125 '^312 = '^312 — 
'^312, we need not to neglect intrinsic piezoelectricity and results are still very simple and 
symmetric: 



h 
h 

''''' ^ ''''' ^ d PoV(^o^33) + eg/Li2 + PoeoMi2' 
^eff ^ , ,h eoPo 

%12 ^312 + ^ p2/(^^^^^) + ,2/,^^^^ + Poeo/rf312' 

IV. DISCUSSION 

Equations ©-(111) express the dependence of sample effective parameters on both material 
and geometric parameters of the layer composite. At first, we can mention that for homo- 
geneous composite, these relations transform to ([S])-©, without the terms containing ^312. 
Intrinsic piezoelectricity is neglected in equations (l2])-(j4]). Second, we can find the effective 
parameters of the ferroic sample with thick metal electrodes and without dielectric passive 
layers, taking the limit of infinite £33^ in (I2l)-(jlj). Third, the denominator on the right in 
discussed relations is the same and we can easily express the ratio of effective parameters. 

Now, we can explore the limits Cq — ^ (pure ferroelectric) or Pq ^ (pure ferroelastic). 
In the first case we get 

^33 = 43^ t/d, (8) 

(1) (2) , 

eff _ -51212-^1212 ^ fQ\ 
^212 - (2) ^ (1) W 
-='1212"' ' -'1212"' 

In equation ([H]), we find the simple result that domainwalls motion is nownot damped 
by mechanical interaction with passive layers, as reported in [gI. In Eq. we can easily 
recognize the pure intrinsic compliance of the composite, without any extrinsic contribution. 
In the second case we get 

e^'^ e^^^ t 

Coo too L 



_eff _ -^33 -^33 /-|r,\ 

33 ~ (2) , , (1), ■ y^^) 

•51212 ~ -5l212 V*^- (l-'-) 
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Now we have got an inverse situation: in equation fllOp . there is pure intrinsic permit- 
tivity, without any extrinsic contribution, while equation ffTTl) is a simple formula that can 



be implicitly found in Naturally, rfgfg = in both limits. We should mention, that 
simple limit — > in equations (IHl) and ( |TT|) does not work, because condition ([T]) for the 
thick sample approximation is violated. In the general case, Pq ^ 0, cq ^ 0, the effective 
permittivity and compliance take the value in the corresponding intervals (fTUl) to ([H]), or ([H]) 
to (HI]). 

Let us now turn to the homogeneous case of equations (jS])-©! where intrinsic piezo- 
electricity is included. First, we can see that intrinsic and extrinsic contributions are now 
explicitly separated. We denote the extrinsic contributions A^o^^, As, Ad. The interrelation 
of the contributions is now rather simple: 

Aeoe : As: Ad = : el : eoPo- (12) 

From equation (12), we get the simple condition 

= 1. (13) 

The authors of Ref. [3] have measured all three effective parameters on the same 
RhH2P0i sample in the temperature interval of 90 to 140 K under the critical point, re- 
sulting in an experimental value of 1,05 for the ratio in equation (fT3|) . They have come to 
the similar result on the base of much more elementary considerations. In |8| , we presented 
discussion of other facts concerning the agreement of theory and that experiment. 

Under the critical temperature of the ferroelectric phase transition, Pq and eo often un- 
dergo substantial changes with changing temperature. Nevertheless, their ratio is almost 
constant in a wide temperature interval. This can be explained by simple physical consider- 
ations. The crystals of KDP family undergo the first order phase transition, which is close 
to the second order one, because discontinuous change of Pq and Cq is relatively small. From 
the Landau- Ginzburg theory that is in a good agreement with experiment in this case, we 
get that both Pq and eo are proportional to (Tc — T)^/^ in the region of temperature T under 
the critical one T^. As a result, the ratios in equation ( |T2l) should remain also constant in 
this temperature interval. This is fully confirmed by the above-mentioned measurements in 
Ref. 13 1 . The situation here is very similar to that discussed in connection with the so- 

n 

called Pippard-Janovec relations (see the excellent theoretical work of Janovec (1966) |14|). 



It is obvious that the results, expressed by equations (fT^ . (|T^ . are general and independent 
on the model. On the other side, results ©-([Tj) are partially characteristic for our model 
because of the presence of geometric factors a and h and the material factors, characteristic 
of the passive layers. 
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